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Abstract 

Background: In order to sustain the gains achieved by current malaria control strategies, robust surveillance 
systems that monitor dynamics of vectors and their roles in malaria transmission over time are essential. This 
longitudinal study demonstrates the trends in malaria vector dynamics and their relative contribution to malaria 
transmission in hyperendemic transmission settings in Tanzania. 

Methods: The study was conducted in two villages within the Kilombero Valley, in rural Tanzania for five 
consecutive years (2008-2012). Seventy-two houses were selected per village and each house was sampled for 
mosquitoes monthly using a CDC light trap. Collected mosquitoes were assessed for species identity and sporozoite 
infection status using PCR and ELISA, respectively. Anopheles funestus and Anopheles arabiensis susceptibility to 
insecticides was assessed using WHO guidelines. 

Results: A total of 100,810 malaria vectors were collected, of which 76% were Anopheles gambiae s. I. and 24% 
were An. funestus. Of all An. funestus samples that amplified with PCR (n = 2,737), 97% were An. funestus s.s., 2% 
were Anopheles rivorulum and 1% Anopheles leesoni. Whereas for An. gambiae s.l. (n = 8,1 1 7), 93% were An. arabiensis 
and 7% were Anopheles gambiae s.s. The proportion of An. gambiae s.s. identified by PCR (2,924) declined from 0.2% 
in the year 2008 to undetectable levels in 2012. Malaria transmission intensity significantly decreased from an EIR of 
78.14 infectious bites/person/year in 2008 to 35 ib/p/yr in 201 1 but rebounded to 226 ib/p/yr in 2012 coinciding 
with an increased role of An. funestus in malaria transmission. Insecticide susceptibility tests indicated high levels of 
resistance in An. funestus against deltamethrin (87%), permethrin (65%), lambda cyhalothrin (74%), bendiocarb 
(65%), and DDT (66%). Similarly, An. arabiensis showed insecticide resistance to deltamethrin (64%), permethrin 
(77%) and lambda cyhalothrin (42%) in 2014. 

Conclusion: The results indicate the continuing role of An. arabiensis and the increasing importance of An. funestus 
in malaria transmission, and pyrethroid resistance development in both species. Complementary vector control and 
surveillance tools are needed that target the ecology, behaviour and insecticide resistance management of these 
vector species, in order to preserve the efficacy of LLINs. 
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Background 

Malaria transmission in humans is sustained through 
vector-human interactions [1] and vector control inter- 
ventions, such as long-lasting, insecticidal nets (LLINs), 
aim to break this interaction. Major promotion of LLINs 
in recent years has resulted in average household owner- 
ship rates and usage of LLINs of approximately 42 and 
36%, respectively, in sub-Saharan Africa [2]. In mainland 
Tanzania, a recent report by the Tanzania HIV and Mal- 
aria Indicator Survey (THMIS) indicates that above aver- 
age LLIN ownership and usage (approximately 90 and 
66%, respectively) was associated with improved malaria 
control and overall reduction in malaria prevalence [3]. 

One outcome of LLIN use is that, by limiting availabil- 
ity of human hosts [4-6], vector species composition in 
any given area can change considerably after a long 
period of LLIN use. Anopheles gambiae sensu stricto, 
Anopheles arabiensis and Anopheles funestus are the pri- 
mary malaria vectors in sub-Saharan Africa [7,8], often 
occurring sympatrically [9]. Anopheles gambiae s.s. is 
often regarded as the most important vector species 
across Africa [9-11] and, because of its almost entirely 
anthropophagic and endophilic behaviour, it is the spe- 
cies that has been targeted most effectively by LLINs. 

However, in some locations, populations of An. gam- 
biae s.s. have developed insecticide resistance and it con- 
tinues to be the dominant vector [12,13]. In other 
locations, An. gambiae s.s. populations have crashed and 
the relative importance of the remaining vector species 
has shifted, with An. arabiensis becoming the major 
malaria vector [4,5,11]. 

Since single populations of An. arabiensis can exhibit 
a range of behaviours, biting and resting indoors as well 
as outdoors and feeding on both humans and animals, 
interventions that optimally target indoor resting and 
biting vectors often impact far less on this species 
[11,14-17]. The primary vector of the An. funestus com- 
plex, An. funestus is also a very anthropophilic and 
endophilic mosquito and it too can be a highly efficient 
malaria vector [10,18,19]. 

Kilombero Valley in southern Tanzania has been sub- 
ject to a large number of studies on malaria epidemi- 
ology, dating back many years, with malaria parasite 
prevalence rates of up to 70% and an entomological in- 
oculation rate (EIR) of 300 infectious bites per person 
per year (ib/p/yr) being recorded in the 1990s, the 
period before the introduction of bed nets [20]. Follow- 
ing the scaling up of untreated nets in the early 2000s 
[21] and insecticide- treated bed nets (ITNs) and LLINs 
from 2004 to 2011 [22-24], a continuous decline in mal- 
aria vector numbers and malaria transmission has been 
seen [3,12]. Although the populations of An. gambiae s.s. 
are significantly dwindling in southern and other parts 
of Tanzania [12], the remaining populations of An. 



arabiensis and An. funestus appears to have shifted their 
blood-feeding periodicity to optimize their chances to 
obtain blood meal from their preferred hosts even in the 
time of low LLIN coverage [25] . It is however suggested 
that prolonged, widespread use of LLINs is likely to 
favour outdoor and early biting, either as an expression 
of the mosquito s innate phenotypic plasticity or possibly 
as a heritable, selectable trait that might be expected to 
increase in frequency [25]. 

The malaria vector populations in this area are subject 
to ongoing rigorous monitoring and herein five years of 
data to the end of 2012 are reported; describing changes 
in vector species composition and relative abundance, 
insecticide susceptibility and their contribution to mal- 
aria transmission following the years of widespread LLIN 
use since first introduced in 2004. 

Methods 

Study site 

The study was carried out in Namwawala (8.154425°S 
and 36.393005°E) and Idete (8.098190°S and 36.510350° 
E) villages (Figure 1) located in the flood plain of the 
Kilombero River (8.1°S and 36.6°E) in southeastern 
Tanzania. The epidemiology of malaria transmission and 
associated vector species composition within these villages 
has been documented over many years [21,25]. Both vil- 
lages experience an annual rainy season (January-May) 
and the main crops are rice and maize. Both villages are 
similar in size (Namwawala = 844 and Idete = 804) and ap- 
proximately 92% of community members sleep under an 
ITN or LLIN [12]. 

Study design 

This study was conducted over five years between Janu- 
ary 2008 and December 2012. A total of 72 houses from 
each village were randomly selected from the Ifakara 
Health Institute (IHI) Demographic Surveillance System 
household list [26]. All selected houses were geolocated 
using a handheld GPS (eTrex, Vista, Garmin, USA). 
Mosquitoes were sampled in every house each month 
during 2008, 2011 and 2012 and for six months from 
January to June in 2009 and July to December in 2010. 

Mosquito sampling and processing 

Mosquitoes were sampled using miniature Center for 
Disease Control (CDC) light traps (model 512, USA). 
One CDC light trap was used overnight per house, 
placed 1-1.5 m from the fan above the ground close to 
the foot end of an occupied bed, and left to run for 
12 hours (19.00-07.00) [27,28]. For every participating 
house, one LLIN (Olyset, A to Z Textiles Mills, Arusha, 
Tanzania) was provided to protect the bed occupant 
where the CDC trap was set. The following morning, 
CDC light traps were collected and mosquitoes killed 
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using chloroform, and identified in the field using a mor- 
phological key [18]. Female mosquitoes were classified as 
being unfed, partially fed, fully fed or gravid. Subsamples 
of five mosquitoes from each trap for An. arabiensis and 
An. funestus species were individually stored inside a tube 
containing cotton wool and silica gel beneath for further 
individual molecular species identification using polymer- 
ase chain reaction (PCR) assay for the An. gambiae com- 
plex [29] and An. funestus group [30] and sporozoite 
infection status using enzyme-linked immunosorbent 
assay (ELISA) [31] in the laboratory (species identification 
for the An. funestus group did not begin until 2009). 

All the sorting information and laboratory analysis re- 
sults were recorded using designated data collection 
forms for entomological studies (Kiware et al, unpub- 
lished). In addition, variations in malaria transmission by 
different vector species over time were assessed and 
compared using the annual EIR calculated by biting rate 
(total collections/trap nights/year) and the proportion of 
females infected with sporozoites [32]. Monthly average 
rainfall data for 2008-2011 were obtained from the Kilo- 
mbero Valley Teak Company (approximately 15 km 
from Idete village), and data for year 2012 data were ob- 
tained using rain gauges installed in Namwawala village. 

Insecticide susceptibility tests 

Following significant increase in An. funestus population 
in 2012, despite extensive usage of LLIN in the study area, 
it was unclear whether this was due to its reduced suscep- 
tibility to the insecticides used in LLINs. The tests were 
conducted using WHO standard procedures and test 
kits for adult female mosquitoes of An. arabiensis and 
An. funestus [33] in Namwawala villages from January to 



June 2013. As the confirmatory process, the biossays were 
repeated in June 2014 for both species. 

Five classes of insecticides currently recommended for 
vector control were tested using discriminating concen- 
trations impregnated in pre-prepared test papers as fol- 
lows: deltamethrin (0.05%), permethrin (0.75%), lambda 
cyhalothrin (0.05%), bendiocarb (0.1%), and DDT (4%). 
Unfed female wild An. funestus collected using CDC 
light traps were used for insecticide exposure bioassays, 
as recommended by WHO for this difficult-to-colonize 
species [33]. However, this method is limited by greater 
variation in susceptibility due to unknown age differ- 
ences between test mosquitoes, it is simple to carry in 
the field with minimal infrastructure and test mosqui- 
toes highly representative of the natural population [33] . 

Prior to exposure, morphologically identified mosqui- 
toes were maintained on 10% glucose solution for at least 
five hours prior to testing; whereas, for An. arabiensis, Fl 
female mosquitoes two to three days old (recommended 
age group) were used for bioassays from reared Anopheles 
larvae collected from the breeding habitats in the study 
sites [34,35]. Species identification was carried out after 
bioassays on dead mosquitoes using PCR. 

A total of 100 mosquitoes were exposed per discriminat- 
ing concentration in five replicates of 20 mosquitoes each, 
and compared to a control with same number of mosqui- 
toes per replicate. In an exposure tube, mosquitoes were 
held for a total of one hour in intervals of 10, 15, 20, 30, 
40, 50, and 60 minutes. After the first hour of exposure, 
mosquitoes were transferred to non-insecticide treated, 
clean, holding tubes and observed for a further 20 minutes 
[33]. After 80 minutes (initial 60 min + further 20 min) of 
knockdown monitoring, all mosquitoes were transferred 
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to non-insecticide treated, clean, holding tubes and kept 
for 24 hours and provided with 10% glucose solution, after 
which mortality was monitored and recorded. All these 
procedures were performed in the field under average am- 
bient temperatures of 26 ± 2°C and a relative humidity of 
78 ± 3% in both bioassay rounds. Percentage knockdown 
in the observed mosquitoes was recorded immediately for 
each time interval, and mosquito mortality in each bio- 
assay was expressed as the proportion of dead mosquitoes 
to total exposed, for each tested insecticide. Execution and 
interpretation followed recently updated WHO test proce- 
dures for insecticide resistance monitoring in malaria vec- 
tor mosquitoes [33]. 

Statistical analyses 

Only data pertaining to An. gambiae si and An. funestus 
were analysed, using SPSS version 20 (SPSS Inc, Chicago, 
USA). Data were fitted with generalized linear models 
(GLMs) using a negative binomial distribution with log- 
link function, and relative rates (RR) with 95% confidence 
intervals calculated to estimate yearly mean mosquito 
catches, relative to the reference year. Species {An. gambiae 
si and An. funestus) were treated as predictors and total 
number of mosquitoes as a dependent variable; the statis- 
tical differences in dependent variables was evaluated as a 
function of villages (Idete and Namwawala), seasons (wet 
and dry) and years (2008-2012). 

Insecticide susceptibility test biossay data were consid- 
ered for each diagnostic concentration and year of testing. 
Mortality was calculated as the percentage of mosquitoes 
dead post 24 hours' exposure to insecticide, and the re- 
sults were assessed according to WHO testing procedure 
for insecticide resistance monitoring in malaria vectors 
[33]. Mortality rates between 98 and 100% indicate full 
susceptibility, 90-97% is suggestive of resistance and re- 
quires further investigation, and mortality rates less than 
90% confirm the existence of resistance. 

Ethical considerations 

The study approval was granted by the Ifakara Health 
Institutional Review Board (IHRDC/IRB/No.A-32) and 
the National Institute of Medical Research (NIMR/HQ/ 
R.8a/Vol. IX/764). On first visiting each house, the bene- 
fits and possible risks associated with the study were ex- 
plained to the house occupants and informed consent to 
proceed was requested. After consenting, the head of the 
house was asked to sign two copies of the informed con- 
sent forms, (retained by the head of the house and the 
study investigator). 

Results 

Relative abundance of malaria vector species 

During the five consecutive years of sampling with CDC 
light traps in sentinel houses, a total of 100,810 malaria 



vectors were collected of which 76% were Anopheles 
gambiae sensu lato and 24% were An. funestus. In each 
of the first four years (2008-2011), the proportion of 
An. gambiae si was significantly higher than An. funestus 
in both study villages (p <0.0001): proportions in total 
catches in Namwawala were 94% (40,028) and 6% (2,398), 
and in Idete were 87% (24,869) and 13% (3,730) for 
An. gambiae si and An. funestus, respectively. However, 
in 2012, the proportion of total catch of An. funestus was 
significantly higher than An. gambiae si. in both villages: 
42% (6,622) and 58% (8,953) in Namwawala and 35% 
(4,479) and 65% (8,447) in Idete for An. gambiae si and 
An. funestus, respectively, (RR (95% CI) = 1.35 (1.23-1.49), 
p <0.0001). 

A total of 8,117 An. gambiae si. were successfully 
identified by PCR and comprised 93% An. arabiensis (n 
= 7,549) and 7% An. gambiae s.s. (n = 568). The relative 
proportions of the species were similar in Idete {An. ara- 
biensis 96% (n = 3,610), An. gambiae s.s. 4% (n = 151) 
and in Namwawala 90% (n = 3,900) An. arabiensis, 10% 
(n = 456) An. gambiae s.s.. However, the relative propor- 
tion between the two sibling species was changing over 
time, with significant decrease of An. gambiae s.s. from 
14% (409/2,924 ) in year 2008 to disappearance 0% (0/ 
1,362 ) in year 2012, compared to An. arabiensis increas- 
ing from 86% in 2008 to 100% in 2012 (Table 1). 

Of the 2,737 An. funestus samples that were identified 
by PCR, 97% were An. funestus s.s. (n = 2,655), 2% were 
Anopheles rivorulum (n = 55) and 1% Anopheles leesoni 
(n = 27). The species composition of An. funestus in 
Idete was 98% (n = 1,554) An. funestus s.s, 1.5% (n = 23) 
An. rivorulum and 0.4% (n = 6) An. leesoni. In Namwa- 
wala it was 98% (n = 1,133) An. funestus s.s., 0.6% (n = 7) 
An. rivorulum and 1.2% (n = 14) An. leesoni, (Table 1). 

Seasonal variation in vector abundance 

During the study, the period from January to May was 
categorized as the wet season, receiving an average 
(+SD) of rainfall of 281 + 178 mm/month, and June- 
December as the dry season, with an average of rainfall 
of 24 + 66 mm/month (Figure 2). The abundance of both 
An. gambiae si and An. funestus peaked in the wet 
season in both villages. The mean number (+SD) of 
An. gambiae si caught per trap per night during the 
wet season was 19 + 48 and 32 + 110, whereas in the dry 
season it decreased to 0.86 + 5.7 and 1.1 + 5.8 at Idete 
and Namwawala, respectively. Furthermore, An. gambaie 
s.s. was only present in the wet season in the first three 
years (2008-2009/10) before its disappearance in 2011/ 
12, compared to its sibling species An. arabiensis, which 
was found to exist in both season, similar to An. funestus 
s.s., a dominating member of An. funestus group. 

The mean number of An. funestus per trap per night 
in the wet and dry season of the first four years of study 
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Table 1 Malaria vector composition, sporozoite 
prevalence (S), biting rate (B) and entomological 
inoculation (EIR) for Anopheles gambiae s.s., Anopheles 
arabiensis and Anopheles funestus and their overall 
estimated yearly contribution to malaria transmission 
from year 2008-2012 in the study area 



Species 


2008 


2009/10 


2011 


2012 


An. gambiae complex sibling species proportion 






An. gambiae s.s. 


0.14 


0.15 


0.002 


0 


An. arabiensis 


0.86 


0.85 


0.998 


1 


No. of PCR amplifications 


2,924 


1,307 


2,542 


1,362 


An. funestus group sibling species proportion 






An. funestus s.s. 




0.887 


0.956 


1 


An. rivulorum 




0.013 


0.021 


0 


An. leesoni 




0 


0.023 


0 


An. parensis 




0 


0.001 


0 


No. of PCR amplifications 




330 


880 


1,527 


Sporozoite prevalence (S; c 


Yo) 








An. gambiae s.s. 


1.18 


0.04 


0 


0 


An. arabiensis 


0.16 


0.36 


0.07 


1.47 


An. funestus 


1.71 


0 


0.43 


2.20 


Biting rate (B; b/p/n) 










An. gambiae s.s. 


8.52 


6.05 


0.04 


0 


An. arabiensis 


52.37 


35.51 


59.74 


20.70 


An. funestus 


1.74 


12.84 


10.09 


14.31 


Entomological Inoculation Rate (EIR; ib/p/y) 






An. gambiae s.s. 


36.70 


1.61 


0 


0 


An. arabiensis 


30.58 


55.51 


15.17 


110.90 


An. funestus 


10.86 


0 


15.58 


115.10 


Total 


78.14 


57.12 


31.05 


226.0 



Note: Sporozoite prevalence = Number of positive sporozoite mosquitoes/total 
tested; Biting rate = Total collections/trap nights/calibration factor, 0.3 for An. 
gambiae complex, and 0.68 for An. funestus [21]; EIR = S x B x 365. 



(2008-2011) was consistently similar in both villages. In 
the wet season, the mean catches (+SD) were 1.23 + 4.7 
in 2008, 2.15 + 7.5 in 2009/10, 0.64 + 1.9 in 2011 com- 
pared to 1.15 + 5.2, 0.77 + 4.3 and 1.62 + 5.52 of the re- 
spective years in the dry season. In 2012, the mean catch 
of An. funestus, both in wet and dry seasons, was ap- 
proximately six times significantly higher than in the 
previous years (p <0.0001): 11.75 + 45.8 and 8.3 + 25.6 of 
wet and dry season, respectively. 

Malaria transmission 

A total of 10,138 individual mosquitoes (530 An. gambiae 
s.s., 7,130 An. arabiensis and 2,478 An. funestus s.s.) were 
screened for Plasmodium falciparum sporozoites of which 
75 were positive (0.74% sporozoite prevalence). Although 
An. gambiae s.s. was the major malaria vector with a 
sporozoite prevalence of 1.18% in 2008, its dominance 



decreased with time to zero in 2011 and 2012, following 
its control to undetectable levels. Conversely, the import- 
ance of An. arabiensis and An. funestus s.s was increasing 
with time from a sporozoite prevalence of 0.16% in 2008 
to 1.47% in 2012 for An. arabiensis, and from 1.71% in 
2008 to 2.2% in 2012 for An. funestus s.s. 

Similarly, the EIR of An. gambiae s.s. decreased dras- 
tically from 30.70 ib/p/yr in 2008 to 0 ib/p/yr in 2012, 
whereas those of An. arabiensis increased approximately 
four times from 30.58 in 2008 to 110.9 in 2012 and that 
of An. funestus s.s. increased 11 times from 10.86 in 
2008 to 115.10 in 2012. 

Overall, the level of malaria transmission in the study 
villages markedly decreased with time from an EIR of 
78.14 ib/p/yr in 2008 to 31.05 ib/p/yr in 2011 but over- 
whelmingly increased to 226 ib/p/yr in 2012, approxi- 
mately seven times more than in the previous year 
(Table 1). 

Anopheles arabiensis and Anopheles funestus insecticide 
susceptibility tests 

In the WHO bioassay testing, as the results indicated 
(Figure 3), An. funestus was fully susceptible to delta- 
methrin (100% mortality) with reduced susceptibility to 
permethrin (93%), and lambda cyhalothrin (91%) and 
confirmed resistance to DDT (86%) in year 2013. In 
2014, An. funestus was resistant to permethrin (65%), 
lambda cyhalothrin (74%), bendiocarb (65%), and even 
to deltamethrin (87%) to which it was fully susceptible 
in 2013. Mortality in control tubes was 4% in both test- 
ing rounds. All tested mosquitoes were amplified as 
An. funestus, using PCR. 

In year 2013, An. arabiensis was fully susceptible to 
bendiocarb (100% mortality) and deltamethrin (98.3%), 
reduced susceptibility against DDT (97%), and con- 
firmed resistance to permethrin (83.3%) and lambda 
cyhalothrin (78%), with a control mortality of 0% across 
all test concentrations. Similar levels of resistance were 
maintained across tested diagnostic concentrations in 
year 2014, whereby the mosquitoes were fully suscep- 
tible to bendiocarb (98% mortality) and resistant to 
deltamethrin (64%), permethrin (77%), and lambda cyha- 
lothrin (42%), with a control mortality of 0% across all 
test concentrations. 

Discussion 

This study provides substantial information on malaria 
vector dynamics and their contribution to malaria trans- 
mission in rural southern Tanzania over a five year 
period. Consistent with other studies, which have docu- 
mented a shift in malaria vector composition and a 
change in malaria transmission dynamics seemingly as a 
result of extensive use of LLINs [4,5,36], this study re- 
ports a steady decrease to undetectable levels of An. 
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Figure 2 Monthly average rainfall in the Kilombero Valley. (A) estimated using CDC monthly biting rates, adjusted by dividing by species-specific 
relative efficiency of 0.3 and 0.68 for An. gambiae s.l. (B) and An. funestus (C), respectively [21], in Idete and Namwawala villages over time. 



gambiae s.s. with steady increase in the proportion of its 
sibling species An. arabiensis and a surge in the abun- 
dance of An. funestus s.s. in year 2012. 

Anopheles gambiae s.s. prefers to feed and rest inside 
houses. This makes it more vulnerable to insecticides 
applied to nets (LLINs) and walls (indoor residual spray- 
ing (IRS)) while An. arabiensis, with its opportunistic 
feeding behaviour both on humans and animals [12,32] 
and its potential to rest outside human dwellings, make 
it less affected by LLINs. Although, lack of outdoor mos- 
quito collections was a major limitation of this study in 
explaining the shift in biting periodicity and outdoor 
biting, it has been recently documented elsewhere that 
An. arabiensis and An. funestus [25,37] display a behav- 
ioural avoidance to contact LLINs by feeding outdoors 
in early part of the evening which might increase its 
chance to survive current interventions. 



A significant increase in An. funestus abundance and 
EIR in 2012 is demonstrated. This shift poses great con- 
cern in malaria control efforts due to its efficiency in 
transmitting malaria. Historically the control of An. 
funestus s.s. was successful through extensive IRS, taking 
advantage of its highly anthropophagic and endophilic 
behaviour, using dieldrin in Pare, Taveta, northern 
Tanzania [38,39] Malindi on the coast of Kenya, using 
DDT [40] as well as in South Africa [41]. This is partly 
because they spend a longer time on insecticide-treated 
materials [42]. However, the vector eventually resurged 
six years later due to a lack of IRS programme continu- 
ity and consolidation [40,43]. A similar scenario was ex- 
pected in this particular region, where usage of LLINs is 
high [3,44]. 

The steady increase in An. funestus population density, 
despite extensive usage of LLINs in the study area, may 
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Figure 3 Results of WHO bioassay test for insecticide susceptibility status of wild female Anopheles funestus (white bars) and Anopheles 
arabiensis (grey bars) from the study sites in the Kilombero Valley, Tanzania, in January 2013 and June 2014. The graph shows percentage 
24 hours mortality rate after a one-hour exposure to the WHO diagnostic doses of insecticide. The minimum sample size for these assays was 1 00. 



be due to its reduced susceptibility to the insecticides 
used in LLINs. Recent findings from western Kenya have 
demonstrated similar phenomenon of resurging An. 
funestus populations, chiefly being due to resistance de- 
velopment to the pyrethroids used in LLINs [45]. 

Preliminary findings from this study demonstrated 
high resistance of An. funestus and An. arabiensis to py- 
rethroids, deltamethrin, lambda cyhalothrin and per- 
methrin, used in Olyset LLINs, distributed in the study 
area in June 2011 [3]. Overall, there was great variation 
of the resistance status between 2013 and 2014 in both 
species tested; however, the variation was surprisingly 
huge in An. funestus than An. arabiensis, which might 
be due to inconsistency in unknown age of the used wild 
mosquito females [33]. 

Due to the absence of organochlorine insecticide DDT 
and carbamate insecticide bendiocarb deployment for 
malaria vector control in the study area, the source of 
resistance in mosquitoes to these insecticides remains un- 
known. Although not tested in this particular study, pyr- 
ethroid (DDT and pyrethroid) carbamate cross-resistance 
was considered to be a probable cause of An. funestus re- 
sistance to DDT and bendiocarb, respectively, which has 
been proved to exist in malaria vectors elsewhere [46,47] . 
In addition, the continuous and illegitimate use of DDT as 
a pesticide in agriculture in the region might have con- 
taminated malaria vector breeding habitats and caused 
physiological resistance in mosquitoes [48]. 



Pyrethroid resistance in both species has been docu- 
mented in multiple countries and regions of East Africa 
[45,49,50], southern Africa [51-54] and West Africa 
[55-57]. Further detailed studies are urgently required to 
establish current vector control operational impacts 
associated with this level of resistance. These findings 
suggest an increased contribution of these vectors to 
malaria transmission and hence great threat to the fu- 
ture use of LLINs in controlling these vectors. 

The other probable cause for the observed increase in 
An. funestus population in this study area, which re- 
quires further investigation, might be a shift of An. 
funestus to outdoor and early evening and daytime biting 
behaviours, which increase their chances to survive and 
reproduce by feeding on unprotected humans, as re- 
cently documented An. funestus behaviours in Benin 
[58] and Senegal [59], West Africa. 

In this study, both An. funestus and An. gambiae s.l. 
vector abundance varied with season. Increases in An. 
gambiae s.l. densities are facilitated by a wide range of 
ephemeral, sunlit, breeding habitats, such as hoof prints, 
rice puddles and ground depressions created during the 
rainy season [18,60]. The temporary nature of these hab- 
itats tends to reduce predation rate but also allows quick 
development of the juvenile stages, which results in 
An. gambiae s.l. dominating during the rainy season 
[18]. On the contrary, An. funestus prefer vegetated 
semi-permanent and permanent breeding habitats, such 



Lwetoijera et al. Malaria Journal 2014, 13:331 
http://www.malariajournal.eom/content/13/1/331 



Page 8 of 10 



as swamps and large ponds [18]. Anopheles funestus 
remained at a reasonable and detectable density across 
the rainy and dry seasons in the study areas and were 
significantly more abundant than An. gambiae si in the 
dry season, probably due to their breeding habitat stabil- 
ity against desiccation [61]. 

Irrespective of seasonal variation in vector abundance, 
An. funestus s.s., An. gambiae s.s. and An. arabiensis 
were all-important malaria vectors in the study area 
[12]. Despite high abundance of An. arabiensis and a 
higher EIR between 2008 and 2010, An. funestus con- 
tributed a relatively higher or equal EIR in 2011 and 
2012. Historically, An. funestus has displayed high sporo- 
zoite prevalence [62] similar to that observed in this 
study and in a recent study conducted in neighbouring 
villages within the valley (Kaindoa et al, unpublished). 
This trend of increase in abundance and high sporozoite 
prevalence of An. funestus has been also observed in 
Asembo district, western Kenya [45] and so appears to 
represent a trend across several regions of East Africa. 

The huge increase in potential malaria transmission in 
2012 (EIR = 226) coincided with an increase in abun- 
dance and sporozoite rates in An. funestus as it did in a 
neighbouring village in the valley (EIR = 467) (Kaindoa 
et al, unpublished). The substantial increase in An. 
funestus and its reduced susceptibility to pyrethroids 
poses a serious threat that needs attention from vector 
control stakeholders. A separate study in West Africa 
also reported a rebound in malaria transmission partly 
being caused by resistance development in An. gambiae 
to pyrethroids [63]. 

A previous study has shown that despite high coverage 
and usage of LLINs, a high proportion of mosquitoes 
still enter houses [64]. Therefore, the increase in An. 
funestus, particularly in the dry season, is likely to ex- 
acerbate the problem. Therefore, new strategies to ad- 
dress resistance and outdoor biting behaviour in the 
early part of the evening as displayed by An. funestus 
and An. arabiensis are required. This can be achieved 
through improving the LLINs; for instance, recent de- 
velopment of nets which can target multiple resistant 
mosquitoes, Olysetf Plus [65], and by targeting vectors 
while outdoors using non-resistant compounds, either 
through larval source management in the dry season via 
autodissemination of insect juvenile hormone, e.g., pyri- 
proxyfen [66,67], or by mosquito sterilization with pyri- 
proxyfen [68], and killing them with toxic sugar-baited 
traps [69], non-chemical electric grid [70] and odour- 
baited traps [71]. 

Conclusion 

This study shows that An. funestus and An. arabiensis 
are important malaria vectors sustaining malaria trans- 
mission, with a substantial increase in An. funestus and 



drastic reduction in An. gambiae s.s. in the year 2012. 
Malaria transmission significantly declined from 2008 
to 2011 and rebounded in 2012 coinciding with an in- 
creased role of An. arabiensis and An. funestus in malaria 
transmission. Although fully susceptible to deltamethrin, 
An. arabiensis and An. funestus were found to be resistant 
and with reduced susceptibility to permethrin pyrethroid 
used for LLINs, respectively. These findings call for com- 
plementary vector control tools, robust vector surveillance 
systems and an insecticide resistance management plan to 
complement and preserve the efficacy of LLINs. 
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